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Abstract 

We explore a mechanism for generating the Yukawa hierarchies by displacing 
the left- and right-handed components of the Standard Model fermions in 
a higher-dimensional space. We find a unique configuration of displacements 
which fits all quark and lepton masses and mixing angles, with a prediction for 
the strange quark mass m^ IS (2 GeV) (1.19) x (V ub V cb /V us ) 1/2 xmf s (m b ) ss 
120 MeV. 
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I. INTRODUCTION 



The parameters of the standard model contain a mystery: why do the Yukawa couplings 
vary over so many orders of magnitude? The up quark Yukawa coupling is of order unity, 
which we would naively expect, but the other couplings are distributed hierarchically over 
several orders of magnitude below this. Why are most of the standard model fermions so 
light? 

Traditional explanations of this Yukawa hierarchy invoke new symmetries, but in |lj 
Arkani-Hamed and Schmaltz suggested an explanation which does not require imposing 
new symmetries. They suggest that the hierarchies present in the Yukawa couplings of the 
standard model can be explained as a result of the slight displacement of the standard model 
field wavefunctions inside a four-dimensional domain wall in a higher-dimensional space. 
The effective four- dimensional Yukawa coupling is a product of the fundamental higher- 
dimensional Yukawa coupling and the overlap of the field wavefunctions. If the wavefunctions 
are highly peaked, a small relative shift between wavefunctions leads to a large suppression 
of the effective Yukawa coupling through the smallness of the overlap of the wavefunctions. 
Even with a single 0(1) higher- dimensional Yukawa coupling, one could produce Yukawa 
hierarchies through the relative displacements of the wavefunction peaks. 

Ref. JIJ also provides a natural, field-theoretic mechanism for producing such displace- 
ments, and demonstrates how it can lead to the appearance of a small Yukawa coupling. 
However, the discussion in |l| is very general, which raises the question of whether a realistic 
model can actually be constructed, or whether there are hidden constraints and relations 
among the masses and mixing angles of the standard model that preclude their generation 
from a set of displaced wavefunctions. 

In this paper we seek to address this question. We find that while it is easy to arrange 
two- (or higher) dimensional arrays of wavefunctions which reproduce the parameters of the 
standard model, it is difficult (but not impossible!) to do so in a minimal one-dimensional 
version; we find essentially one possible configuration of wavefunctions. This configuration 
has a definite prediction, which may be confirmed or falsified by more accurate determina- 
tions of the strange quark mass and V u t,. For alternative extra-dimensional approaches to 
flavor see 

In section [H], we summarize the mechanism of and in section [TTI| we summarize 
the outline of a model suggested in that paper. In section [iy| we demonstrate the mech- 
anism for the simple case of the lepton sector, and in section ^ we outline the chal- 
lenges present when trying to model the quark sector. In section [VT] we describe our 
efforts to find a set of locations for the quark wavefunctions, and in section |VII| we ana- 
lyze the (essentially) unique set, which leads to the prediction for the strange quark mass 
m s m (1.19) x rribCVubVcb/Vus) 1 / 2 « 120 MeV. We briefly conclude in section |VTT|. Appendix 
[A| contains a derivation of the form of the zero mode for a minimal five-dimensional spinor 
coupled to a scalar with a domain wall profile expectation value, and appendix O contains 
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the experimentally allowed masses and mixings we used, as well as the renormalization 
multipliers used to evaluate those parameters at a common scale. 

II. PRODUCING HIERARCHIES WITH EXPONENTIALLY SMALL OVERLAPS 

In this section, we summarize the mechanism of |T[] for producing hierarchies between 
mass matrix elements. We start by considering a theory in which the standard model 
fermions are trapped in a four-dimensional membrane in a five-dimensional space. Although 
non-perturbative effects in string/M theory can be used to achieve such localizations, one 
can also study models in which the localization occurs in the context of an effective field 
theory. 

One concrete field-theory method of achieving such a localization is to couple a massless 
five-dimensional fermion f to a five-dimensional scalar field $ through a Yukawa coupling 
term ~ / d 5 x§^ty in the Lagrangian, and then to include an effective potential for $ with 
two or more isolated vacua, such as a potential ~ ($ 2 — v 2 ) 2 , so that it attains a position- 
dependent vacuum expectation value with a domain-wall profile in the extra dimension, 
interpolating between regions with different choices of vacuum (see figure |l|). 

The fermion then acquires a zero-mode localized near the zero-crossing of $. Heuristi- 
cally, one may think of the Yukawa coupling to the position-dependent scalar expectation 
value as giving the fermion a position-dependent mass which greatly suppresses field fluctu- 
ations far from the domain wall at £5 = while leaving fluctuations at the origin massless. 
Although a single minimal spinor in five-dimensional space decomposes into a left-right 
mirror pair of chiral fermions in four- dimensional language, only one chirality zero-mode is 
trapped on the wall. 

In section ^ we derive the zero mode solution in the approximation in which the scalar 
field profile is a linear function of the extra dimension <3> = 2/z 2 xs (which is valid for points 
close to the center of the domain wall) and find that zero mode solutions for the fermions 
have Gaussian profiles in the extra dimension centered on the zero of 

V(x) = A e~^ xs2 x i[>(x Q , x u x 2 , x 3 ), (2.1) 

where ip is a canonically normalized massless left-handed four- dimensional fermion field, 
A = ij}I 2 {ix /2)~ l l A is a normalization constant, and \x is related to the slope of the scalar 
field profile by /i 2 = <9 5 ($)/2. Higher modes, which come in mirror pairs, have minimum 
energy « \i. 
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FIGURES 




FIG. 1. A chiral fermion trapped on a domain wall. 

In [[!]], it was pointed out that adding a five-dimensional fermion mass term / d 5 xM^ty 
to the above system acts like an effective shift of the scalar field $ — ► $ + M as far as the 
fermion is concerned, which (in the approximation in which the scalar field profile is linear) 
leads to an effective shift of the location where the scalar has a zero expectation value, i.e., 
a shift in the effective location of the domain wall as seen by the fermion (see figure |2|) so 
that it ends up localized around £5 = — M/2/i 2 . If different fermions have different masses 
Mi, they each end up localized around different locations £^5 = — Mj/2/i 2 , each fermion 
being centered on a different "effective" zero of the same scalar field. 




FIG. 2. A mass term shifts the location of the domain wall. 
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We assume a generic five-dimensional Yukawa coupling between two fermions and a Higgs 
scalar field H, which we take for simplicity to have a constant profile in the extra dimension. 
This leads, upon substitution of the zero mode solution from equation ( |2.1| ) and integration 
over the extra dimension, to a four dimensional coupling modulated by the mutual overlap 
of the fermion wavefunctions: 



d 5 x 


\TL k i 


d 5 x 


Vl k 


d 4 x 


(Jdx* 


d 4 x 


( e -^(i 


d 4 x 


A hiptfa 



C Y uka W a = I d b x VL K H^ 2 (2.2) 

\= hj ( A e-^>-^ 2 ^) (A e-^-^fo) (2.3) 
i e-^^-^A e -» 2 (xt>-t*A hfafo (2.4) 
] \) hfaifo (2.5) 

(2.6) 

where we have inserted a factor of the square root of the domain wall width \fh to render 
the coupling constant k dimensionless. The effective four-dimensional Yukawa coupling: 

A = e-W^K, (2.7) 

depends on the mutual overlap of the zero modes of the fermions through the relative 
displacement of their peaks \£i — If the fermions are more than a few Gaussian widths 
from each other then an order-one coefficient for the five-dimensional Yukawa coupling 
k results in an exponentially small four- dimensional Yukawa coupling A. Setting the Higgs 
field to its expectation value and generalizing to many generations gives a mass matrix 
element: 

my = ((h) k) e-^Hi-ij) 2 

= pe -W&-W\ (2.8) 

where p = (h) k. In this paper we seek to obtain the hierarchical pattern of masses and 
mixings in the standard model as a result of a particular configuration of wavefunction 
locations with relative separations of a few Gaussian widths yT 1 . 

It is useful to invert equation Q2.8| ) to find the relative displacements associated with a 
mass matrix: 



= rf2 1og(^/p) = p(mj)- (2.9) 

To produce a particular mass matrix m^-, we need to find a set of wavefunction locations 
which have relative displacements \£{ — £j\ related to the mass matrix by the above formula. 
The function p~ l ^—2 log(my/ p) appears often in expressions to follow; we refer to this 
function as p(rriij). 
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III. THE MODEL 



We now summarize the model outlined in Q. By the standard model gauge symmetry, 
all fields in a given gauge multiplet have the same effective five dimensional fermion mass 
Mi and thus have wavefunctions peaked at the same location in the extra dimension. This 
means that there are fifteen wavefunction locations qi,Ui,di,li and e% for the left-handed 
standard model fields Qi, Uf, Df, Li and Ef (with i = 1, 2, 3). The fermions are all taken to 
have equal couplings to the domain wall scalar $ (meaning their Gaussian wavefunctions all 
have equal widths), and so as not to add extra complication, we take the Higgs field H to 
have an expectation value which is constant throughout the domain wall, and we take the 
effective five dimensional Yukawa couplings between the fermions and the Higgs to have a 
single magnitude. Although these restrictions could be relaxed, the structure already present 
proves sufficient to generate the standard model masses and mixings. We defer a discussion 
of the CP violating complex phase of the CKM matrix to later work, and for the present 
time study only real Yukawa matrices. 

The analysis of Arkani-Hamed and Schmaltz constrains the scales present in the model. 
For the description to make sense, we need the wall thickness L to be larger than the 
Gaussian width fj," 1 , which in turn should be larger than the length scale of the ultraviolet 
cutoff 1/M*. For the four dimensional effective top Yukawa coupling At to be perturbative 
at M* we must further constrain NXf/Wn 2 < 1 where N = (M 3f L)/2n is the number of 
Kaluza-Klein modes below the cutoff. If a field theory description is to apply throughout 
the wall, we must have $(L/2) ~ p?L < M#. To summarize: 

L- 1 < [i < M* (3.1) 
MJL- 1 < 32n 3 /\ 2 t « 1000 (3.2) 
fx/L- 1 < y/MjL- 1 « 30. (3.3) 

Note that these inequalities only constrain the ratios of the mass scales. The overall 
mass scale corresponds to a free parameter of the model. It would be most interesting if this 
overall mass scale were low, so that the extra-dimensional physics would be accessible to 
experiment. Constraints on the size of the extra dimension L come from direct searches for 
Kaluza-Klein excitations and precision electroweak data H|. Constraints on split fermions 
were discussed in However, the most important constraints come from considering 

flavor changing neutral currents mediated by Kaluza-Klein gauge bosons || , and lead to the 
bound 

L- 1 > 100 TeV . (3.4) 

In equation ( |2.8| ) the mass matrix element is bounded from above by p, with the maximal 
value achieved when the fermion wavefunctions are localized around the same point. This 
means that to generate the top mass from a single matrix element the condition p > m t is 
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required. On the other hand, naturalness and perturbativity of couplings lead us to want 
p as small as possible. For concreteness, we choose p to be 1.5 m t . Our model works for 
any p > m t , with larger p values requiring greater suppressions from the overlaps, and thus 



is independent of p, since the wavefunction positions are always chosen to reproduce the 
same mass matrices. 

Before trying to reproduce the masses of the standard model, we first run the masses of 
the standard model to a common scale m t using the scaling factors found in the appendix 
(the running of the CKM elements to that scale is small compared to the experimental 
uncertainty in their magnitudes). 



We begin our efforts to reproduce the standard model masses and mixings by finding a 
set of wavefunction positions for the Lj and Ef fields which reproduce the observed charged 
lepton masses. This is easy, since we are only trying to fit three parameters and we may 
adjust five parameters (six locations, minus one freedom to shift the locations by an overall 



We have three constraints on the locations, with the only other constraints being that the 
relative displacements indicated by (far) be great enough to lead to negligible off-diagonal 
matrix elements. In fact, it turns out that little care is necessary to insure that the off- 
diagonal elements are negligible; one need only insure that, for a given L field, the closest 
E c field is roughly one width closer than the second-closest E c field. It proves to be a general 
property of this model, due to the exponential in equation ( |2.8| ), that if one wavefunction is 
near two others and the two others are not the same distance away then the overlap with 
the farther one will often be negligible compared to the overlap with the closer one. In other 
words, it is easy to generate elements of the Yukawa matrix that are effectively zero. 

One possible set of locations for the lepton fields is: 



requiring larger separations between wavefunctions. The mass prediction of equation ([7.131) 



IV. LEPTON SECTOR OF THE MODEL 







(4.3) 
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This configuration is depicted graphically in figure [3[ This is by no means the only set 
of wavefunction locations that reproduce the e, /i, and r masses. For instance, the mass 
spectrum is essentially unchanged if one moves the Li further apart while keeping each 
Ef wavefunction the same distance from its partner Lj. Also, one can consider textures 
differing from that of equation (|4.1| ) — since lepton mixing angles are unobservable, all that 
is necessary is that the eigenvalues of the mass matrix equal the lepton masses. Different 
textures would lead to different sets of wavefunction locations. 



L wavefunction 
E wavefunction 

x m u e e 

L ' L 'R R I 
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FIG. 3. The locations of the lepton wavefunctions. 



V. QUARK SECTOR OF THE MODEL 

We now outline the challenges present when trying to model the quark sector. It is more 
difficult than the lepton sector because one has to fit nine observables (six masses and three 
mixing angles) by adjusting eight parameters (the nine locations of the Q iy U?, and D°, 
minus one freedom to shift the locations by an overall displacement). 



In fact, equation (2.8) determines the mass matrices m u and m^. To obtain the physical 
masses and CKM elements from the mass matrices we must perform unitary redefinitions of 
the Q (1=4-1/2)*) Q(i=-i/2)i, U? and D\ fields to a basis where the mass matrices are diagonal. 
The masses can be read off the diagonal entries, and the CKM matrix appears in the weak 
interaction terms as a result of the mismatch between the redefinitions of the isospin-up and 
isospin-down components of the quark doublets. 

A useful method of computing the parameters is to calculate: 



/ml 

m\ | (5.1) 

V m\ 
(m\ \ 

m d mi = U d ml U d (5.2) 

V ml) 

Vckm = UtU d . (5.3) 
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The eigenvalues of the hermitian squares of the up (down) mass matrices give the squares 
of the masses of the up (down) type quarks, and the products of the eigenvectors of the 
hermitian square of the up mass matrix with the eigenvectors of the hermitian square of 
the down mass matrix give the CKM matrix elements. Ignoring the CP violating complex 
phase, it is sufficient to match the magnitudes of the three components of the CKM matrix 
above the diagonal to match the entire matrix. 

However, we seek to do more than simply reproduce the standard model masses and 
mixing angles — we seek to reproduce the hierarchies present in a natural way, as a result of 
the exponential in equation ( |2.<j| ). This excludes certain mass matrix textures from consid- 
eration. 

For example, consider the mass matrices: 

and i'l 'I). (5.4) 




Both yield a spectrum with a hierarchy (in fact, they each yield a spectrum containing one 
massless and one massive particle) but in the first matrix the hierarchy is a result of the 
hierarchy between the last element and all the other elements, which equation ( |2.8| ) easily 
produces, while in the second matrix the hierarchy is a result of a delicate cancellation 
between all the matrix elements, which would be exponentially fine-tuned if these elements 
came from equation (p^8|). Although these fine tunings could possibly be explained by adding 
further symmetries, our philosophy is to seek a set of generic-looking positions which lead to 
small masses and mixings as a result of small numbers coming from the exponential in Q2.8|). 



We only allow small numbers to come from the product of small numbers, but not from a 
cancellation or sum. One might be tempted to allow cancellations or sums from numbers 
which are of the same order-of-magnitude as the quantity they cancel or sum to produce, 
since it does not seem like fine-tuning to have 0(e n )±0(e n ) = 0(e n ), but we avoid doing this 
here because in this model it would usually involve an unexplained coincidence — generically 
any two elements of the mass matrix resulting from equation ( |2.8j ) have vastly different 
orders of magnitude and it takes a strange coincidence or a fine tuning for two elements to 
have an 0(1) ratio. An exception to this is the third generation: the masses of the top and 
bottom approach the upper-bound scale p, and as we approach this scale such coincidences 
become less fine-tunedQ so we allow models which produce a mass parameter ~ m t or ~ mj 
through the sum or difference of two terms of roughly that order of magnitude. We first 
consider the case where m t comes from a single mass matrix element and m?, also comes 
from a (different) single matrix element (we will later need to relax the constraint on m^). 



1 The fine tuning drriij/dij — > as \£i — ► 0, which happens when rriij — > p. For \£i — £j\ larger 
than a few widths fj, , that is, for rriij below the mass scale of the third generation, drriij/dlj 
grows rapidly. 



9 



A simple mass matrix texture appropriate to our philosophy is: 



ra, 



m d 






(m d 








VcKM ' 





m s 







V o 





m b 



(5.5) 



(5.6) 



in which each small entry arises directly from the exponential in equation 

We now try to construct these matrices from a model. Translating the matrices in Eqs. 
( |5.5| ) and ( |5.6| ) into matrices of constrained relative distances using equation ( |2~9"| ) gives: 

(far) 

Vh - " ,1 = l-r' [jar) 3.4840 (far) ] (5.7) 

(far) 
4.4278 



\Qi - dA 



-i 




4.0715 

(far) 



(far) \ 
(far) 
.90052/ 
4.5319 s 
3.9378 
3.0029 



(5.8) 



These constraints can easily be satisfied in two-dimensional extensions of the current model 
involving string defects in place of domain walls (see figure |j). Unfortunately, the configu- 
ration cannot be collapsed down to one dimension without bringing some wavefunctions too 
close to preserve the smallness of the "zero" elements of the texture. 




10 12 



FIG. 4. A possible set of locations for quark wavefunctions in a model with two extra dimen- 
sions. The lines represent constrained distances between wavefunction locations. 
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A similar texture can be found by taking the down mass matrix diagonal and using: 

/ m u m c V cd m t Vtd ' 



m, 



m c m t V t s 
V m t 





/ m u 







~ v CKM 





m c 


: ) 




V 








(5.9) 



and, at the same level of approximation, one can use textures where two of the Vqkm 
elements come from one mass matrix and one comes from the other. The textures can be 
further extended by allowing the small masses to result from products of two small mass 
matrix entries by generalizing the observation that the 2x2 mass matrix: 

(5.10) 




with b,c <C D has masses D and be/ D. For instance, the mass matrix: 

m s V us m h V uh 

m d = I m d /V us m s m b V cb I , (5.11) 




m b 



leads to a down quark mass m d . 



VI. SEARCHING PARAMETER SPACE 

In this section we describe our efforts to find a set of locations for the quark wavef unctions. 
The textures described in section [V] (where each mixing parameter comes from a nonzero 
above- diagonal element in either the up mass matrix or the down mass matrix, and the 
masses come from either a diagonal element or a product of a below-diagonal element and 
an above-diagonal mixing term) were searched analytically with no success. The closest 
match allowed one to fit all the masses and two of the mixings by hand, with all the "zero" 
elements of the texture being of negligible magnitude. Fitting these parameters determined 
all the wavefunction positions, and thus predicted the remaining mixing (parameterized by 
V u b) which was found to be many orders of magnitude too small. We then resorted to 
numerical methods, and also relaxed our assumptions to allow for m t to result from the 
contributions of two or more elements of roughly equal magnitude (and later extended this 
relaxation to nib)- 

We implemented a brute-force scan over the parameter space, which yielded a solution 
with two mass matrix elements of magnitude m b described in section |VII| below. 



VII. ANALYSIS OF SOLUTION 

We now describe the configuration of quark wavefunction locations which was found 
to produce agreement with the observed masses and mixings of the standard model. The 
configuration: 
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-i 



-7.6057\ 
6.9522 ] , 
0.0 



-2.7357X 
10.4362 
0.9012 / 



11.3682 \ 
-3.2250 
3.0511 / 



(7.1) 



was found to produce mass matrices:^ 



m, 



/ 1.7630 
1.0365 x 10~ 15 

5902.8 
m u 
m c 
m t . 



5.1637 x 10~ 66 

576.06 
5.5689 x 10~ 19 



4.8087 x 10- 11 
2.7882 x 10~ 3 
165900 



MeV 



(7.2) 



(7.3) 



and: 



m d 



1.6660 x 10~ 73 

14.510 
2.1526 x 10~ 23 



16.947 
8.0344 x 10~ 18 
1373.2 

m s V us /cos{6) 
m d /V us m b V cb /cos(6) 

m b sm(6) m b cos(9) 



5.4422 x 10- 20 
123.42 
2370.2 



MeV 



(7.4) 



(7.5) 



which, using equations ( |5.1| ) - ( |573| ) and running using the scaling parameters in the appendix, 
gives masses: 



m u = 3.24 MeV 
m c = 1.25 GeV 
m t = 166 GeV 

and CKM matrix elements: 



m d ■ 

m s 

m b 



6.01 MeV 
■ 120 MeV 
4.25 GeV 



(7.6) 



V 



CKM 



/ 0.9755 0.2200 0.0031 
0.2197 0.9748 0.0390 
\0.0116 0.0373 0.9992, 



(7.7) 



which are consistent with experiment. The configuration is depicted graphically in figure |5|. 



2 The element m u n 3) = 5902. 8MeV is negligible because it is in the same row as the much larger 
top mass, and can be removed by a small 1-3 rotation redefining the Uf fields. 
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FIG. 5. The locations of the quark wavefunctions. 



The down mass matrix texture of equation (|7.5|) differs from the example texture of 
equation ( |5.11| ) only by a 2-3 rotation by 9 re-defining the D\ fields, given that a certain 
additional relation among the entries holds. The rotation on the D\ fields transforms the 
down mass matrix from equation ( [7.5|) as: 



m d — > m c 



/ m s V us m s V us tan(0) 

m d /V us m b V ch tan(0) m b V cb j . (7.8) 
V m b 



which agrees with the texture of equation ( |5.11|) if we can make: 

m s V us tan(6 l ) = m b V ub (7.9) 
m b V cb tan(6 l ) = m s . (7-10) 

Since 9 is a free parameter we can always choose it to satisfy one of the two equations, but 
the other will only be satisfied if 

m b V ub m s 

— — = tan (9) = — — . 7.11 

m s V us m b V cb 

As the mass of the strange quark is the most uncertain of the quantities in equation ( [7.1 1[) , 
we re-write this relation as: 

m s = (WK,) 1/2 xm,, (7.12) 

which holds at the common scale m t . Running the masses down to their physical scales 
using the factors in the appendix gives: 

mf(2 GeV) « (1.19) x (V ub V cb /V us ) 1/2 x mf(m t ) « 120 MeV. (7.13) 

The relation of equation (|7.13 ) is consistent with experiment, and represents the "prediction" 
present in this model which allows it to fit nine observables with eight model parameters. 
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There are also seven other configurations giving the same masses and mixings which result 
from the fact that, for all three i, moving the Uf fields to the opposite side of its partner 
quark doublet Qi while keeping their relative distance fixed leaves the diagonal elements of 
m u unchanged while not carrying the Uf so close to the other quark doublets Qj^i that they 
create large off-diagonal elements in m„. As these other configurations produce the same 
mass matrices, they yield the same prediction for the strange mass of equation ( |7.13| ). 

No other such distance preserving re-arrangements to the configuration can be made 
without bringing otherwise distant fields too close and thus creating large unwanted mass 
matrix elements which destroy the predictions above. The constrained distances between 
the wavefunction locations are illustrated in figure []. 
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FIG. 6. The constrained distances between the quark wavefunctions. 



VIII. CONCLUSION 

It is pleasing to find that Arkani-Hamed and Schmaltz's solution to the Yukawa hierarchy 
problem can be made concrete, and that the masses and mixings of the standard model, 
while highly constraining, do not have any hidden relations that prevent their generation 
from a configuration of wavefunction displacements. The minimal model described contains 
the prediction Eq. ( [7.13Q which may be falsified or supported by future experimental data. 

It is a pleasure to thank Nima Arkani-Hamed, Travis Brooks, Lance Dixon and Michael 
Peskin for useful discussions and suggestions. This work was supported by the Department 
of Energy under contract DE-AC03-76SF00515. 

APPENDIX A: DERIVATION OF THE ZERO MODE 

In this appendix we demonstrate that a minimal five-dimensional spinor coupled to a 
scalar field with a generic domain-wall expectation value profile in the fifth dimension leads 
to a single chiral fermion zero mode localized at the zero crossing of the scalar field. We 
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show that if the scalar profile is linear, the zero mode has a Gaussian profile, with the higher 
modes having the spectrum and profile associated with harmonic oscillator wave-functions. 
The five- dimensional action: 



(Al) 



S = J d 5 x ^[ij A d A + $(ar 5 )]tf, 
can be decomposed into four-dimensional terms: 

S = [d 4 x [ dx 5 Uli-fdp + z 7 5 <9 5 + $(x 5 )]*. 



In chiral spinor notation, this is: 



S 



d 4 x / dx 5 [V>1 iPr] 





-<9 5 + $(x 5 )~ 






+d 5 + $(x 5 ) 






_tpR_ 



(A2) 



(A3) 



We can split the spinors into a 4-dimensional spinor part ip multiplied by an extra- 
dimensional profile (with no spin indices) £: 



and integrate over x$ to get: 

S = 



E 



^L(i){x,,)C L {x 5 ) 



(A4) 



(A5) 



where the m; and C/'' are chosen such that: 



(_a 5 + $(x 5 ))Cf (x 5 ) = m 4 C L (x 5 ) 

(+9 5 + $(x 5 ))C i ^5) =m 4 Cf(x 5 ), 



and we have defined: 



From equations 



with solutions: 



and 



*0*W(i) 

we find zero modes: 



;-9 5 + *(x B ))tf(x B ) = 
(+9 5 + $(x 5 ))Co L (^) = 0, 



Cf(a; 5 ) = <e + r 5 *fc)^ 
C L (o; 5 ) = ^e-r^)^, 



(A6) 
(A7) 



(A8) 



(A9) 
(A10) 



(All) 
(A12) 
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where the A are normalization constants chosen to make / dx 5 (*( = 1. In the approximation 
$(£5) = 2fi 2 xs this is: 



C B (x 5 ) = Ae + ^ (A13) 
C i (x 5 ) = Ae-^i, (A14) 

so the left-handed solution Co( x 5) is a Gaussian and the right-handed solution Cq(xs) is non- 
normalizable — only the left-handed solution is physical. By inspection of equations (|A11 ) 



and ( A12 ) we see that this is a generic feature of any function $(^5) that goes from a finite 



negative asymptotic value to a finite positive asymptotic value with one zero crossing. If the 
Yukawa coupling were multiplied by —1 or if the scalar profile went from positive to negative, 
the right-handed solution would have been the physical solution and the left-handed solution 
would have been non-normalizable. 

If the scalar expectation value dipped back through zero, as in a configuration of two 
domain walls coming from a kink/antikink profile, there would be a mirror zero-energy chiral 
fermion stuck on the other wall. In the limit that the second domain wall (and its mirror 
fermion) are infinitely far away, we recover the case discussed above. 

To examine the massive spectrum, we operate on equation (|A6|) with (+d$ + $(xs)) and 
equation ( |A7|) with (— d 5 + $(x 5 )) and use $(x 5 ) = 2/i 2 x 5 to obtain: 



+ ( 2/ u 2 x 5 ) 2 + 2/, 2 )Cf (x 5 ) = m 2 Cf (x 5 ) (A15) 
(-«9 2 + (2/i 2 * 5 ) 2 - 2/i 2 )C L (x 5 ) = m 2 C L (x 5 ). (A16) 

(A17) 

Recognizing these as harmonic oscillator Schrodinger equations for a particle with mass 
h 2 /2, characteristic angular frequency ujq = 4/i 2 /fr and ground state energy ±2/z 2 , we can 
immediately write down the solution: 

m 2 n = Afi 2 n {n = 0,1,2,3,...) (A18) 

cfe) = c? +1 (* 5 ) = (^rcfe) - 4t ( ~ 95 t Vx5 ) n fe). (A19) 

In this framework we can see the chirality of the zero mode by noting that the ground state 
energy shift term ±2fi 2 of the Schrodinger equation shifts the eigenvalue tower from the 
usual ~ 1/2,3/2,5/2... down to ~ 0,1,2... for the left handed modes and up to ~ 1,2,3... 
for the right handed modes. 

APPENDIX B: STANDARD MODEL PARAMETERS 

For the experimentally allowed values of the standard model parameters, we used [0: 
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m u = 1.5 to 5 MeV 
m d = 3 to 9 MeV 
m s = 60 to 170 MeV 



m u /m d = 0.20 to 0.70 
(m u + m d )/2 = 2 to 6 MeV 

V ;/ = 34 to 51 

m d - m u 

m c = 1100 to 1400 MeV 
m b = 4100 to 4400 MeV 
m t = 166000 ± 5000 MeV 



V us = 0.217 to 0.224 
V ub = 0.0018 to 0.0045 
Va = 0.036 to 0.042 



m e 



0.5110 MeV 
105.7Mey 
1777 MeV 



(Bl) 



where the quark masses are in the MS renormalization scheme. The up, down, and strange 
masses are evaluated at a scale of 2 GeV, while the other quark masses are evaluated at a 
scale equal to their MS mass. The lepton masses are pole masses. 

To run the quark masses to the common scale m t , we used the three- loop QCD and 
one-loop QED scaling factors & 



Vu 


Vd 


Vs 


Vc 


Vb 


Vt 


1.84 


1.84 


1.84 


2.17 


1.55 


1.00 



(B2) 



where Vi = m i{ m i) / m i{ m t) for % = c,b,t and t]i = m,i(2 GeV)/mi(m t ) for % = u,d,s. The 
lepton pole masses were related to their MS masses evaluated at m t using the relation 

,2 



rn 



pole 



m 



MS 



M 1 + 



a 



71 



1 + ?log(^) 



to give scaling factors 



Ve 




Vt 


1.05 


1.03 


1.02 



(B3) 



(B4) 



where r\i 



mf ole (mi ) / mf IS (m t ) 
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